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Abstract 

The cosmological singularities of the Bianchi I universe are analyzed in the setting of loop 
geometry underlying the loop quantum cosmology. We solve the Hamiltonian constraint of the 
theory and find the Lie algebra of elementary observables. Physical compound observables are 
defined in terms of elementary ones. Modification of classical theory by holonomy around a loop 
removes the singularities. Ho"wever, our model has a free parameter that cannot be determined 
■within our method. Testing the model by the data of observational cosmology may be possible 
after quantization of our modified classical theory. 
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Typeset by REVI^ 



1 



I. INTRODUCTION 



The Bianchi I universe is of primary importance as it underlies, to some extent, the 
Behnskii-Khalatnikov-Lifshitz (BKL) scenario HI, [s], 0, @], which is beheved to model 
the Universe in the vicinity of the cosmological singularity. It has been studied recently 
[sl, Isl, [lo|, 11, 12] within standard loop quantum cosmology (LQC). 



The standard LQC [H, M means basically the Dirac method of quantization, which 
begins with quantization of the kinematical phase space followed by imposition of constraints 
of the gravitational system in the form of operators at the quantum level. Finding kernels of 
these operators helps to define the physical Hilbert space. In the nonstandard LQC [H, 1^ 
one first solves all the constraints at the classical level to identify the physical phase space. 
Next, one identifies the algebra of elementary observables (in the physical phase space) and 
finds its representation. Then, compound observables are expressed in terms of elementary 
ones and quantized. The final goal is finding spectra of compound observables which are 
used to examine the nature of the big-bounce phase in the evolution of the Universe. 

This paper is devoted to the classical dynamics of the Bianchi I model with massless 
scalar field modified by the loop geometry, described in the framework of the nonstandard 
LQC. The next paper will address the problem of the quantum dynamics [l7j . 

In Sec. II we define the modified classical Hamiltonian. Section III concerns the final 
choice of canonical variables. The classical dynamics is solved in Sec. IV. An algebra of 
elementary observables is the subject of Sec. V. Physical observables, that may be confronted 
with the cosmological data, are defined in Sec. VI. We conclude in the last section. Appendix 
A presents derivation of the symplectic form on the constraint hypersurface. In Appendix 
B we derive the algebra of elementary observables corresponding to the Bianchi I model 
without the loop geometry modifications. 



II. HAMILTONIAN 

The gravitational part of the classical Hamiltonian, Hg, in general relativity is a linear 
combination of the first-class constraints, and reads 

Hg := / d^'x^N'Ci + N^'Ca + NC), (1) 

where S is the spacelike part of spacetime M x S, (iV*, iV°, N) denote Lagrange multipliers, 
{Ci,Ca,C) are the Gauss, diffeomorphism and scalar constraint functions. In our notation 
(a, 6 = 1,2,3) are spatial and {i,j,k = 1,2,3) are internal SU{2) indices. The constraints 
must satisfy a specific algebra. 

The Bianchi I model with massless scalar field is described by the metric: 

3 

ds"^ = -N^ dt^ + ^ a]{t) dxl (2) 

i=l 

where 

k * 3 3 

aiir) = ai{0) (-] , dr = N dt, ^ k; = 1 = ^ kf + kj, (3) 

i=i i=i 
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and where describes matter field density (k^ = corresponds to the Kasner modep. For 
clear exposition of the singularity aspects of the Bianchi I model we recommend 0, l^- 

Having fixed local gauge and diffeomorphism freedom we can rewrite the gravitational 
part of the classical Hamiltonian, for the Bianchi I model with massless scalar field, in the 
form [3] 

Hg = -7"' / d^x Ne-'e,,kE^^E'"'F^, , (4) 
Jv 

where 7 is the Barbero-Immirzi parameter, V C S is an elementary cell, S is spacelike 
hypersurface, N denotes the lapse function, Sijk is the alternating tensor, E°- is a densitized 
vector field, e := a/ | det E \ , and where is the curvature of an SU{2) connection A^. 

The resolution of the singularity, obtained within LQC, is based on rewriting the curvature 
F^ij in terms of holonomies around loops. The curvature may be determined by making 
use of the formula 

= -2 lim Tr ( ~ ] "ui, (5) 

where 

is the holonomy of the gravitational connection around the square loop Oij, considered over 
a face of the elementary cell, each of whose sides has length ^jLj (and Vo := L1L2L3) 
with respect to the flat fiducial metric °qab '■= ^ij°'-^a°'-^i] the fiducial triad °el and cotriad 
°uj'^ satisfy °ujl°e'j = 6y, ArOij denotes the area of the square; and = /y y/°qd^x is the 
fiducial volume of V. In what follows, to simplify further discussion, we make the assumption 
Li = 1 = L2 = -L3 (so Vo = 1), which is natural in the case of choosing the Bianchi I model 
with T^-topology instead of M^. 

The holonomy in the fundamental, j = 1/2, representation of SU(2) reads 

h'f'^ = cos(/iiQ/2) I + 2 sin(/i,Q/2) r,, (7) 

where = — ■icrj/2 (cxj are the Pauli spin matrices). The connection and the density 
weighted triad (which occurs in f lTU]l ) are determined by the conjugate variables c and p: 

A\ = d Et = p. (8) 

where: 

Ci = 7di, \pi\ = aj Ofc (9) 
Making use of (jl]), ([H]) and the so-called Thiemann identity 

c- ^-1 T^aj T^bk _ Sgn(piP2P3) ^ o^abco, ,krj.^( (Tn\ 



k 



leads to Hg in the form 



where 



Hg= lim /7(a'ia'2M3)^ (11) 



(12) 
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and where V = aia2 03 is the volume of the elementary cell V. 

The total Hamiltonian for Bianchi I universe with a massless scalar field, 



reads 



H = Hg + Ha, 



0, 



(13) 



where Hg is defined by (111 I) . The Hamiltonian of the scalar field is known to be: = 
N p'^\p\~^ /2, where and Prp are the elementary variables satisfying {4>,p^} = 1. The 
relation H ^ defines the physical phase space of considered gravitational system with 
constraints. 

Making use of ([7]) we calculate (I12p and get the modified total Hamiltonian Hg'^^ corre- 
sponding to (1131) in the form 



1 



Sgn(piP2P3) 



sin(ci/ii) sin(c2Ai2) At3 sgn(p3 



\PlP2\ 



+ cyclic 



+ 



where 



\Pi 



A, 



(14) 
(15) 



and where A is a regularization parameter. Here we wish to emphasize that (1141) presents a 
modified classical Hamiltonian. It includes no quantum physics! 

In the gauge = P2 Psl the Hamiltonian modified by loop geometry reads 



^(A) 



87rG'72A2 

The Poisson bracket is defined to be 



|piP2|^^^ sin(ciyUi) sin(c2/i2) + cychc 



d- d- d- d- 



+ 



d- d- 



d- d- 



'dp^ dp^d(f)' 



^dckdpk dpkdck 

k=l 

where (ci, C2, C3,pi,p2,P3, 0,P(/>) are canonical variables. The dynamics of ^ reads 

i ■= {^,H^^^}, ^ e {Ci,C2,C3,Pi,P2,P3,0,P<^}- 



(16) 



(17) 



(18) 



The dynamics in the physical phase space, ^^yg, is defined by solutions to (|T8l) satisfying 
the condition H^^"^ ^ 0. The solutions of (fT8|) ignoring the constraint H^'^^ ^ are in the 
kinematical phase space, J-'k^l- 



III. NEW CANONICAL VARIABLES 

We use the following canonical variables 

Ci 



P^■.-- 



Vi ■= |pi|^''^sgn(pi 



\Pi 



where i = 1,2, 3. They satisfy the algebra 

{Pi,Vj} = UnGjSij, 



(19) 



(20) 
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where the Poisson bracket reads 

3 



{■,-} = 12vrG7 



d- d- d- d- 



dPkdvk dvkd(3k 



k=l >- 

The Hamiltonian in the variables ffT9l) turns out to be 



d- d- d- d- 

+ TT 



rr(x) pI 1 /sin(A/3i) sin(A/32) , sin(A/3i) sin(A/33) 

= TT - n ^ 9 ^1^2 H 715 ^^l^'S 



2 87rG72 V A2 ' " ' A2 

sin ( A/32 ) sin (A/33) 
+ V2V3 

where A parametrizes the holonomy of connection modifying the Bianchi I mode 



IV. DYNAMICS 

A. Equations of motion 

The Hamilton equations of motion read^ 

P. = -18vrG^^(0, + 0,), 



Vi = 187rG't;iCos(A/3i) (Oj + Ok) 

= P^, 

P<f> = 0, 

H^^^ ^ 0, 



where 

^ _ Vi sm{Xf3i 



B. Solution to equations of motion 

Insertion of into (123|) gives 

_ _tan(AA) 

A Vi 

which leads to 

sin(A/3i) 

Vi ; = const 

A 

Therefore, Oj are constants of motion. 



^ where i,j,k = 1,2,3 and i ^ j ^ k 



Makiner use of ([25]), ([21]) and cos(AA) = - sm(AA)2 erives 



^V^^-i^TlG^X Oi)2 

Integration of fl3T]) leads to 



ISvrG 



(Oj + Ok) 



(31) 



In 



ISvrG 



(Oj + Ok) (0- 



Thus we have 



2|t;,| =exp ( (Oj + Ok) 



X exp 



- 0°) ) + (127rG7AOi)2 x 
ISttG 



(32) 



(33) 



(Oj + Ok) 



°)p 



which may be rewritten as 



Vi = 127rG7A O; cosh ( (Oj + Ok) 



6°) -ln|l27rG7AOi| 



(34) 



It results from the above solutions that for a nonzero value of A there is no Big Bang 
type singularity (for any value of 0). The Big Bang is replaced by the Big Bounce. In [7| 
one considers the so-called planar collapse, but we do not consider this issue here as we are 
mainly concerned with an initial type singularity. 

Using flMl) it is not difficult to get 



sin(A/5j 



1 



cosh ' ^^'"^ 



p,j> 



(Oj + Ok) 



^) -ln|l27rG7AOi 



(35) 



Alternatively, one may solve the equation of motion to get f3. From (l25l) and ( l23l) we obtain 



A dpi 



-18nG 



The integration of ([36]) gives 



In 



tan 



(Oj + O, 



(Oj + Ok) 



(36) 



+ const. 



(37) 



Removing the cosmological singularities does not complete our task. In what follows we 
consider an algebra of elementary observables and phy sical compound observables as they 
define the background of the nonstandard LQC (l5l.1l6|. 
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V. OBSERVABLES 



A function F defined on the phase space is a Dirac observable if it is a solution to the 
equation 

{F,/7(^)}^0. (38) 

An explicite form of (l38l) is given by 



127rG'7 J2 



1=1 



OF dH^^^ dF dH^^^ \ OF 



d(3i dvi dvi d(3i ) 



+ 



0, 



(39) 



which due to (22) reads 



i=l 



dF sin(AA) dF 
..cos(AA)^ 



dF 

■ (Oj + Ok) + = 0. 



(40) 



A. Kinematical observables 

One may easily verify that Oi satisfy ( HOj) . Instead of solving (140|) one may use the 
constants that occur in flM|) and fl5S]) . This way we get 



Ai =ln 



tan (^) 



+ ISvrG 



(Oj + Ok) 



(41) 



The observables (HT!) are called kinematical as they are not required to satisfy the con- 
straint (ETD. 



B. Dynamical observables 

An explicite form of the constraint fl27|) in terms of Oi is given by 

sgn(p^) = 6v^ VO1O2 + O1O3 + O2O3. 



(42) 



It results from (22), (27) and (28) that O1O2 + O1O3 + O2O3 > so (42) is well defined. 
Thus, the dynamical observables, Af^"^, corresponding to ( 14T1) read 



Af^" = In 



tan (M) 



3v/^sgn(p^)(0j + 0k) 
VO1O2 + O1O3 + O2O3 



(43) 



C. Algebra of observables 

One may verify that Af^'^ satisfy the following Lie algebra 

{Oi,Oj} = 0, 
{Af^^O,} = 5.,, 
{Af^",Af°} = 0. 



(44) 
(45) 
(46) 



In the physical phase space the Poisson brackets are found to be (see, Appendix A) 



^ f d- d- d- d- \ 
-Idyn := g - dO^^) ' ^^^^ 

and the algebra reads 

{Oi,Oj}dyn = 0, (48) 

{Af^^Oj}dyn = 5,,, (49) 

{Af^°,Aj^^^}dyn = 0. (50) 

VI. COMPOUND OBSERVABLES 

In what follows we consider the physical observables which characterize the singularity 
aspects of the Bianchi I model. It is helpful to rewrite (H2l) and flM|) in the form 



pI = SGTrG (O1O2 + O1O3 + O2O3) , (51) 
^ 12.G,X I 0.1 cosh pv^sgn(,,)(O, + O.)0 / A^ _ 

' ' " V VO1O2 + 01O3 + 02O3 V2; w ^ ' 

The so-called directional energy density |3| is defined to be 

p,{\<P):=^,. (53) 

The bounce in the z-th direction occurs when pi approaches its maximum 0], which 
happens at the minimum of Vi [p^f, is a constant of motion). One may easily verify that in 
the case when all three directions coincide, which corresponds to the Friedmann-Robertson- 
Walker (FRW) model, these densities turn into the energy density of the fiat FRW with 



massless scalar field 15 



It is clear that Vi takes minimum for cosh(-) = 1 so we have 



12.G,A0.. = (54) 



Rewriting Oi and p^ in terms of ki and fc^ [3] 



2 



Oi = - kiir, p^ = VSvrG k^K, (55) 



where is a constant, leads to 

We can determine p'^°'^ if we know A. However, A is a free parameter of the formalism. Thus, 
finding the critical energy densities of matter corresponding to the cosmic singularities of 
the Bianchi I model is an open problem. 
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One may apply (l56i) to the Planck scale. Substituting X = Ipi gives 



p-- ^ 0, 35 (^-^ J ppi, (57) 

which demonstrates that p™-"^' may fit the Planck scale depending on the ratio k^/ki. 

Another important physical observable is the volume of the Universe. From the definitions 
(19]) and ([T9D we get 

V = 010203 = \viV2V3\^^^. (58) 

It results from (152|) . (!55l) and ([3]) that the volume is bounded from below. 



VII. CONCLUSIONS 

The modification of the classical Hamiltonian by using loop geometry turns classical 
singularities of the Bianchi I model into Big Bounces, similarly as in the case of the initial 



singularity of the FRW type models [15 



Our approach is quite different from the so-called effective or polymerization method (see, 
e.g. [3]), where the replacement (3 sin(A/5)/A in the Hamiltonian finishes the procedure of 
quantization. In our method this replacement has been done entirely at the classical level. 
Quantization consists in finding a self-adjoint representation of observables on the physical 



phase space and an examination of the spectra of these observables |15l . |16 . 

The elementary observables constitute a complete set of constants of motion on the 
constraint surface. They are used to parametrize the physical phase space and are "building 
blocks" for the compound observables like the directional energy density and the volume 
operator. So they have deep physical meaning. Their role becomes even more important at 
the quantum level as they enable finding quantum operators corresponding to the classical 
compound observables [16]. 

The modification of classical theory by holonomy of connection around a loop does not 
change the Lie algebra of elementary observables. The algebras of modified fj48|) - fl50|) and 
nonmodified flB9l) - flBll|) observables are, to some extent, isomorphic. It is a valuable feature 
of the modification procedure. 

Our quantum Bianchi I model has a free parameter A that cannot be determined within 



the model. This parameter, similarly to the case of the FRW [15|, is expected to be fixed 
by the data of observational cosmology (see, e.g. [l8|). 

The algebra of elementary observables is defined on the physical phase space. The com- 
pound observables are thus defined on the physical phase space too. Thus, their properties 
may be confronted with the data of observational cosmology, i.e. real world. 

The maximum of the energy densities may fit the Planck scale which seems to be surpris- 
ing as the equation used concerns only the classical level. But a similar situation occurred 



in the case of the FRW model [15|] and the explanation has come from the examination of 



the quantum energy density operator. Its spectrum turned out to coincide with the classical 



counterpart |16 



The present paper gives a background for quantization of the Bianchi I model [17[. The 
quantization is required despite the fact that the singularity problem is resolved already 
at the classical level due to the modifications based on the loop geometry via holonomy of 
connection. It is so because the spectra of the quantum observables may be used to get a 



9 



link with the data of observational cosmology, similarly as in the case of the FRW model 
Another reason is that the Big Bounces may occur at any energy density (being 



parametrized by a free parameter) so it is necessary to quantize the model. 

Our paper shows that the nonstandard LQC method, worked out for the case of the FRW 
model [l5|, may be applied to the Bianchi I model. It seems to be applicable to the isotropic 
models (e.g. Lemaitre) as well. 

Acknowledgments 

We are grateful to Przemyslaw Malkiewicz and Wojciech Zaj§czkowski for helpful discus- 
sions. WP thanks Haiyun Huang for interesting correspondence. 

APPENDIX A: SYMPLECTIC FORM 

The symplectic form on the physical phase space, fi, may be obtained from the symplectic 
form on the kinematical phase space, u;, by taking into account the constraint (H2l) . 
The symplectic form corresponding to (12T|) reads 

1 ^ 

uj = V dpi A dvi + d(pA dp^. (Al) 

i/7rG7 ^ — ' 

i=l 



Makin use of the constraint ( 1421) gives 



V f ^ + 1^ dv}\ = t 3V^^g^(^^)(Q. + Q0 ,0. (A2) 



Thus, we have 



^ J. 3v/^sgn(^,)(0, + 0,) = V ^ ^0 A c^O.. (A3) 



VO1O2 + O1O3 + O2O 



3 



i=l 



On the other hand (j28]l leads to 



which gives 



Thus, we have 



dvi = ^^^.^J"^ (iOi - cot(A/5i) (A4) 
sm(A/^i) 



dPi A (it;^ = "''^^,^3'^ dPi A c?Oi. (A5) 
sm(A/^i) 



1 3 ^ A 

-— ^ VrfAAdi;, = V ■ ,.^. c;AAdO,. (A6) 
IzttGj ^ ^—^ sm(Apj) 



It results from ( ]A3[) and (lASp that we have 



3 / \ (9A'^-^'^ 
^ Vsm(AA) 90 / 
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Now, let us rewrite dAf^^ as follows 
Next, we make summation 



Eventually, we calculate the wedge product with dOi 

E ^A?- A = g (-^ + ^ d^. j A (AlO) 
Finally, the physical symplectic form Q reads 

3 

= ^ dA^y'' A dOi. (All) 

Thus, the physical phase space may be parametrized by the variables Af^"^ and Oi, and the 
corresponding Poisson bracket is given by fH7j) . 

APPENDIX B: NONMODIFIED CASE 

To find an algebra of elementary observables, for this case, we introduce the definitions 
Oi : = lim Oi = (Bl) 

Ai-. = lim Ai = In [A[ + ISttG ^ ^^'^ 0. (B2) 



A-»0 



In the limit A the Hamiltonian constraint (|5T1) turns into "unmodified" constraint 



sgn(p<^) = 6v^ ^/O^O^TO^O^TO^s. (B3) 
It results from (lBl])-(lB3]) that 



.dyn ^ In I /J I , Sgn(p^)(C'j + Cfc) (/ 



(B4) 



One may verify that 

{0„0,} = 0, (B5) 
= (B6) 
= 0, (B7) 
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and 



Direct calculations lead to 



n = J2 dAt""" A dO^. (B8) 



i=l 



{Oi,0,}dyn = 0, (B9) 

{^f^^O,}dyn = (BIO) 
{^f^°,4^^}d,n = 0, (BU) 



which coincides with the algebraic structure of the algebra of the modified observables ME 

m. 
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